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We use a simple a transverse radial boost scenario coupled to PYTHIA events to illustrate the 
impact radial flow may have of two- and three-particle correlation functions measured in heavy-ion 
collisions. We show that modest radial velocities can impart strong modifications to the correlation 
functions, some of which may be interpreted as same side ridge and away side structure that can 
mimic conical emission. 
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I. INTRODUCTION 

Observations of Au + Au collisions at the relativistic heavy ion collider (RHIC), and comparisons with results from 
p+p and d+Au interactions have led to the notion that strongly interacting quark gluon plasma is formed in Au + Au 
collisions. This conclusion is in part based on the observation of strong elliptic flow [l|, dramatic reduction of the high 
transverse momentum (pt) particle production (relative to expectations from p + p interactions), and suppression of 
the high transverse momentum two-particle back-to-back correlations ■ Observations of Au + Au collisions have also 
lead to discoveries of new, rather unexpected features, which await satisfactory interpretation. Among these, we note 
the observation of a flattening, or formation of a dip on the away-side in two particle azimuthal angle correlations 0], 
and the discovery of a strong ridge-like structure on the near-side of two particle correlations measured in relative 
pseudo-rapidity vs azimuthal angle space [3| ■ The away side dip and ridge are illustrated qualitatively as cartoons in 
Figure [TJ 





FIG. 1: (Color online) Cartoon illustrations showing (left) the away-side dip structure observed after elliptic flow subtraction 
in A(t> correlations, and (right) the ridge structure observed on the near-side in At; vs. A(f) correlations. 



The observation of a dip at 180° in flow subtracted two-particle azimuthal correlations measured in Au + Au 
collisions was reported by both STAR [1, Q, and PHENIX [J. Stocker argued in [1] the propagation of jets through 
matter, and associated energy losses may generate shock fronts and deflected jets. Shock fronts and wake fields 
calculations by authors of f9l.1l0l. [ll| suggest the observed away-side structures might indeed result from a Mach cone 
effect. Heinz et al. however have argued, on the basis of 2-D hydro calculations, the jet energy loss process is unlikely 
to result in massive away-side structures [l2| . Other works suggest the dip- like structure might result from large angle 
gluon radiation [l^, [3l , Cerenkov gluon radiation [l^, [l^ , or medium drag effects [l^l ■ Note also studies based on the 
AMPT model which suggest that parton transport may only partly explain the observed two-particle correlations [3 ■ 
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There are thus many suggested explanations for the observed away-side structures found in two-particle correlations. 
It is conceivable that discrimination of these different production mechanisms may not possible with two-particle 
correlations, it might however be achieved with three particle correlations. Three particle correlation measurement 
techniques and related issues have been discussed by one the authors and others [13, [21I. [2p. Three-particle 
correlation studies aimed at resolving the above production mechanisms are works in progress j2lll23l [2^ . 

The near side ridge-like structure was first reported by the STAR collaboration [2^, and has been subject of 
extensive studies The ridge yield for a given associate particle pt range is essentially independent of the trigger pt- 
Additionally, particles in the ridge exhibit pt spectra slopes that are independent of the trigger particle momentum. 
This suggests that while the ridge accompanies jet-like particles, it may not be strictly caused by jets, but "associated" 
with jets. This point of view is perhaps reinforced by the observation that the near side zt di-hadron fragmentation 
function, subtracted for ridge contribution, is independent of the trigger momentum [J. 

The near-side ridge has also been subject to numerous theoretical investigations. Interpretations of the ridge include 
induced radiation coupled to longitudinal flow [26| , turbulent color fields [27j , anisotropic plasma [28j , a combination 
of jet quenching and strong radial flow [2^, and recombination of locally thermal enhanced partons produced by jet 
energy loss [sO]- A recent study based on the AMPT event generator [ll| indicate significant string melting though 
it produces a ridge like structure cannot account for the longitudinal extent of the ridge reported by STAR. 

Additionally, it is important to realize that at relatively large pt, the away-side jet is found to reappear. Indeed, 
STAR reported already at QM05 that away side high pt particles accompany trigger particles with pt > 8 GeV/c [sif . 
The away side yield measured in central Au + Au is relatively weaker than that found in d -I- Au, but the azimuthal 
profile of the high pt away- side is otherwise unmodified relative to d -I- Au collisions [3^ . 

It is thus clear that the structure of jets as measured through two-particle correlations is significantly modified in 
central Au + Au collisions. The cause of the modifications and the emergence of new features such as the away-side 
dip and the near-side ridge are however as of yet unclear. While many authors have argued the newly found structures 
arise from jet energy loss and coupling of the jets with gluons of the medium, we will here argue that these may be 
associated with radial flow effects. Evidence for elliptical flow is now firmly established, and it is widely interpreted as 
resulting from low viscosity hydrodynamics [s^ . Radial flow must also be present given the strong elliptical flow. In 
fact, particle spectra, as well as HBT correlations are well described by Blast-wave parametcrizations which include 
large transverse flow velocities (see for example (s^ . 
correlations have already been discussed by Voloshin 
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In this work, we study the effects radial flow can impart on two- and three-particle correlations associated with jets, 
and their underlying events. We neglect string melting, jet quenching, and diffusion processes, and use in-vacuum 
p-l-p events, simulated with the PYTHIA event generator, to study the effect of radial flow on produced jets, and the 
2-/3-particle correlations they induce. We boost particles produced by PYTHIA radially in the transverse plane, at 
flxed velocities, and in an arbitrary azimuthal direction. We argue that a form of radial flow drives the formation of 
the near-side ridge and away-side dip in nuclear collisions. The elliptic and radial flow of particles at soft scales is 
experimentally well established and well described by hydrodynamics (37l . [ssj . The elliptic flow measured in lepton 
production indicates that charm quarks also flow; see [s^, ISl and references therein. It is therefore not unreasonable 
to explore how flow might affect the jet production. 

To understand how jets can acquire a flow-like radial motion, consider a quark antiquark pair produced in a hard 
scattering in a nuclear collision. The center of mass of the pair can have a net transverse momentum at the time it 
is produced due to initial state parton scattering. In essence, the partons that will ultimately produce the qq pair 
can scatter before the jets are produced.. Such scattering takes place in the production of lepton pairs by the DrcU 
Yan process [4l|, and is well understood on the basis of pQCD [4^. This mechanism explains the observed increase 
of (pt) for lepton pairs produced in pA collisions [i^, and also contributes to the Cronin effect in jets [3]. We would 
expect this contribution to be most pronounced in the pt region where the Cronin effect is important. If initial state 
scattering is the primary source of "jet flow" , then we could expect to see appropriate angular correlations in hard 
lepton pairs in pA collisions. 

We view initial state parton scattering as the least speculative mechanism that may contribute to jet flow. Additional 
radial motion can be acquired before the hard-scattering if radially-directed color flelds are established very early in 
the collision, as considered by Fries, Kapusta and Li [4^ . This contribution need not be limited to the Cronin region. 
The fluid flow of the medium can also add momentum to quark and antiquark pair if the hard scattering that produces 
the qq occurs sufficiently far from the center of the collision volume [46|. For example, flow would push both the 
q and q radially if they are produced back-to-back in a direction perpendicular to the radius. Jet quenching would 
accentuate this effect in the correlation function by suppressing the contribution of pairs in other conflgurations (i.e., 
with either the qorq moving radially inward. If this were the dominant mechanism, then we expect the jet flow effect 
to be correlated with the flow plane. 

The mechanism we discuss here should apply to quark jets of flavors, as well as gluon jets. Electron measurements 
suggest that charm quarks indeed exhibit elliptic flow; see [l^SO] and references therein. This suggest that flow could 
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also contribute to jet modification. 

As a starting point in this work, and to emphasize the potential role of flow on two particle azimuthal correlations, 
we show in Figure [5] results from a di-jet toy model in which the near and away side jets are simulated with Gaussian 
profiles of 0.4 and 0.5 radians widths respectively. Jet particles are simulated as pions with transverse momentum 
slope of 300 MeV/c, and have relative probabilities of 0.7 and 0.3 on the trigger and away sides respectively. Figure 
[2] shows correlation functions obtained with flat radial boost profiles for (3r < 0.1 in black, 0.2 < Pr < 0.3 in red, and 
0.4 < (3r < 0.5 in blue. The (3,. < 0.1 distribution features peak structures on both the near and away side, consistent 
with di-jet events. One finds however the introduction of even modest boost (0.2 < j3r < 0.3) produces a flattening 
of the away side jet peak, whereas stronger velocities result in a wide dip on the away-side as observed in STAR and 
PHENIX data. 

The above toy model is rather primitive, so we turn to PYTHIA to provide more realistic jets, and underlying 
events. PYTHIA offers the advantage of providing well tested descriptions oi p+p events and jets. It includes effects 
of momentum conservation, and resonance decay known to influence the shape and strength of correlation functions. 
A key question is whether, or how, jets couple to radial flow. In this work, we will assume maximum coupling: Jets 
produced by PYTHIA will be used, without quenching or modification, to model the effect of radial flow on the shape 
and strength of two- and three-particle correlations. 




FIG. 2: (Color online) Two-particle azimuthal correlations obtained with a di-jet toy model assuming Gaussian azimuthal 
profile and flat radial boost profiles - described in the text. 

This work is divided as follows. Section 2 presents definitions and examples of two- and three-particle correlation 
functions based on PYTHIA p + p collisions generated at ^/s = 200 GeV. In Section 3, we introduce transverse radial 
flow and show by comparison to reference plots of Section 2 how even modest flow velocities imparted to PYTHIA 
events can qualitatively reproduce the away-side dip and near-side ridge structures observed experimentally. 

II. REFERENCE CORRELATIONS 

We first present reference correlations based on p-l-p collisions at y/s = 200 GeV simulated with PYTHIA, version 
6.4, using default parameters (i^. PYTHIA (versions 6.2 and 6.3) simulations were found to qualitatively reproduce 
inclusive and identified particle spectra measured by STAR in p + p collisions at ^/s = 200 GeV [H, |4§| , as well 
particle yield ratios [H^. Figures [H O and [7] present two-particle correlations in A77 vs. A(f> while figures 01 [SI and 
|8] display three-particle correlations in A(j)i2 vs Ai^is relative azimuthal angles. All correlations arc calculated for 
charged particles produced in the range \r]\ < 1. Two-particle correlations are studied between "high" and "low" 
momenta particles. For three particle correlations, we use one high momentum "trigger" particle, and two lower 
momenta particles in the same range. Specific transverse momentum ranges used in the different correlations are 
discussed next. 

We illustrate that the strength and shape of the correlations depend on the transverse momentum ranges considered. 
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Panel (a) of figures 3, 5, and 7, show two-particle correlations obtained for a jet tag particle in the range 3 < pt < 
20 GeV/c, and "associates" in the ranges 0.2 < pt < 1 GeV/c, 1 < pt < 2 GeV/c, and 2 < pt < 3 GeV/c, respectively. 
The top subpanel of each figure show contour plots of the normahzed two-particle density, l/N2°'^(PN2/dAridA(j), 
where A^|°* is the integral of the correlation function. The middle and bottom subpanels display projections dN2/dAr] 
and dN2/dA(t> obtained for |A0| < 0.7 rad, and IA77I < 1.0, respectively. The vertical axes of the projections are 
labeled A'2 for brevity in the figures. 

The two-particle densities, and their projections, obtained for different pt ranges exhibit similar features but with 
varying strengths and widths. The peak centered at A77 = 0, A(/) = corresponds to the "triggered" or tagged jet. 
The narrow ridge-like structure centered at Acj) = tt which extends over the full A?7 = 2 range covered by the plot 
arises from the fact that in p + p intereactions, jets are produced by partons colliding with typically different fractions 
"x" of the proton momenta. Consequently, while the jets are nearly back-to-back in azimuth (because of momentum 
conservation) , they may have quite different rapidities thereby leading to an elongated away-side structure as seen in 
panel (a) of figures [31 and [T] 




FIG. 3: (Color online) Normalized two-particle density, l/Ni°*d'^N2/dAr]dAcf) obtained witli p + p PYTHIA events; Particle 
1 (jet tag particle) in the range 3 < pt < 20 GeV/c, particles 2 and 3 (associates) in the range 0.2 < pt < 1 GeV/c. Panel 
(a) through (d) obtained with radial boosts /3r=0.0, 0.2, 0.3, 0.4 respectively. Middle and bottom subpanels show projections 
along Arj for \A(t>\ < 0.7 rad, and for \Ari\ < 1.0. See text for details. 

It is also useful to consider the jet correlation structures measured on the basis of three-particle azimuthal corre- 
lations to shed some light on features expected in ongoing search for conical emission in heavy ion collisions. Panel 
(a) of figures SI [HI and [5] display three-particle correlations in the form of 3-cumulants plotted as a function of two 
relative azimuthal angles A^i2 and Aipi^. The middle and bottom subpanels show projections of the 3-cumulant 
along diagonal axes. The middle panels projections are carried along the 3-cumulant main diagonal A(/)i2 -I- Aipi^ with 
a condition \ A(j)i2 — A^iaj < 0.35, and shifted so their origin coincide with the away-side (A0i2 -I- A(/)i3)/2 — tt = 0. 
The bottom panels show projections of the 3-cumulant along the second diagonal A(j)i2 ~ A(j)i3 with the condition 
|(A(/)i2 + A(/)i3)/2 — 7r| < 0.3 chosen to emphasize the shape and strength of the away-side peaks. The 3-cumulant 
used here was introduced in reference [l^, and is defined as 

n A J, \ (12) (3) (13) (2) (23) (1) , „ .-.s 

C3(A0i2, A(^13) P3 - p\ - P2 p\ - 92 'Pi'+2pipipi (1) 
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FIG. 4: (Color online) Normalized three-particle azimuthal cumulant C3(A(^i2, Ac^ia) obtained with p+p PYTHIA events: 
Particle 1 (jet tag particle) in the range 3 < pt < 20 GeV/c, particles 2 and 3 (associates) in the range 0.2 < pt < 1 GeV/c. 
Panel (a) through (d) obtained with radial boosts Pr—O-O, 0.2, 0.3, 0.4 respectively. Middle and bottom subpanels show 
projections along cumulant diagonals. See text for details. 



in terms of the measured 3-particle density (number of triplets/event) = dN^ / dA(j)i2dA(f>i3 and combinatorial 
terms P2''^ p^i^ and pipipi calculated based on measured 2- and 1- particle densities. The 3-cumulant is used, instead 
of a simpler 3-particle density ^3^^^"* = (i^iV3/o?A(/)i2(iA(/)i3, because it enables a measurement of three correlated 
particles by eliminating the combinatorial contributions 

pi'\f^ where p^'^'^ = dN2/dA4>ij and ^ = dNi/d4>k. 
The 3-cumulants shown in panel (a) of figures [H El [SI are obtained with the same pseudo-rapidity, and transverse 
momentum cuts as those used in Figures [31 [5l [T] The high momentum cut defines particle " 1" , while the low 
momentum cut is used to select particles "2" and "3". The 3-cumulants exhibit peak structures expected for two jets 
produced nearly back-to-back in azimuth [l^. The peak centered at A0i2 = 0, A(/)i3 = corresponds to two "low" 
momentum particle being emitted along side the "trigger" high pt particle. The peak centered at A0i2 = tt, A0i3 — it 
corresponds to two "low" momentum particle being emitted nearly back-to-back with the "trigger" high pt particle, 
and thus correspond to an away-side jet. Peaks at A(/)i2 = 7r(0), A(/)i3 = 0(7r) correspond to combinations where 
one of the low momenta particle is along-side the high pt particle and the other back-to-back. One observes from 
panel (a) of figures 01 [51 [5] that the strength of the 3-cumulant decreases progressively for increasing selected associate 
momenta. One additionally finds that peaks are narrower at high momentum owing to the fact that high pt particles 
tend to be better aligned with the direction of the jet they are produced from. Also note that the away side peak 
(A</)i2 = TT, A(j)i3 = n) is wider along the main diagonal than the alternate diagonal. This is usually associated to 
the fact that the incoming partons that produce the jets are not co-linear in the laboratory reference frame, i.e. they 
carry a small but finite transverse momentum: the jets are therefore not exactly back-to-back in azimuth, thereby 
leading to finite elongation of the peak along the main diagonal A0i2 -I- A^ia axis. 
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FIG. 5: (Color online) Normalized two-particle correlations in A77 vs. Ai^ obtained with p+p PYTHIA events: Particle 1 (jet 
tag particle) in the range 3 < pt < 20 GeV/c, particles 2 and 3 (associates) in the ranges 1 < pt < 2 GeV/c. Panel (a) through 
(d) obtained with radial boosts /3r=0.0, 0.2, 0.3, 0.4 respectively. Middle and bottom subpanels show projections along A77 for 
IA0I < 0.7 rad, and A(j> for |Ar;| < 1.0. See text for details. 



III. THE EFFECT OF RADIAL FLOW ON TWO- AND THREE- PARTICLE CORRELATIONS 



We now proceed to study the impact of radial flow on two- and three- particle correlations. We use single p+p 
PYTHIA events radially boosted at fixed velocity in random directions in the transverse plane to simulate the effect 
of radial flow on particle correlations. The use of a single event presents the advantage one does not need to carry 
background subtraction of combinatorial terms that arise in A-l-A collisions. One can picture the phenomenon along 
two different scenarios: first consider two incoming partons are subject to transverse kicks before they collide to 
produce jets. As stated in the introduction, this phenomenon has been invoked successfully to interpret Drell Yan 
data [4l[, and also contributes to the Cronin effect in jets In this scenario, the jet production in the rest frame 
of the parton pair is essentially the same as if the partons were not radially boosted, but the produced hadrons are 
boosted in the lab frame thereby leading to modified correlations as shown in the following. 

In an alternative scenario, one considers incoming partons interact as if in vacuum, but their collision products 
are subject to a transverse boost due to momentum kicks by the dense medium. This assumes jet hadronization 
proceeds quickly within the medium, and sufficient medium fiow builds up to radially push hadrons produced by 
parton fragmentation. One does not expect all hadrons to be subjected to the same kick but overall to a finite average 
kick. 

We assume that in either case, the particle production, locally in the parton-parton rest frame, proceeds similarly 
as in vacuum, so as to preserve hadron correlations. While we do not, in practice, expect the partons to behave quite 
that way, the above scenarios enable a simple exploration of the effect of radially boosted particles produced a by 
reasonably well calibrated common source: PYTHIA. We note however that the boosting of an entire PYTHIA event 
rather than only hadrons from jets results in an artificial enhancement of the near correlations. 

In a nucleus-nucleus collision, we would further model the radial boost velocity as being proportional to some 
small power of the distance, in the transverse plane, between the parton-parton interaction point and the nucleon 
participants center-of-mass position. However, one should also account for quenching of the jets, as well as rescattering 
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FIG. 6: (Color online) Normalized three-particle azimuthal cumulant Cs{A<j}i2, Acpn) obtained with p+p PYTHIA events: 
Particle 1 (jet tag particle) in the range 3 < pt < 20 GeV/c, particles 2 and 3 (associates) in the range 1 < pt < 2 GeV/c. Panel 
(a) through (d) obtained with radial boosts /3r=0.0, 0.2, 0.3, 0.4 respectively. Middle and bottom subpanels show projections 
along cumulant diagonals. See text for details. 



and diffusion effects. Such complications are here neglected by considering correlation induced by one radially boosted 
p+p collision only. Note additionally that in an actual A + A collision, multiple p+p collisions would suffer a radial 
boost in different directions. As these are to first order uncorrclatcd, this would lead to a large flat background in 
the correlations under study. The strength of the correlations above this background shall depend on the number of 
p + p 01 "di-jet" sources. Such effects are also neglected. 

We study the impact of radial boost in finite increments of velocity to understand the modifications imparted to 
the particle correlations. Panels (b), (c), and (d) of figures 4-8 show two- and three- particle correlations obtained 
with radial velocities /3,. =0.2, 0.3 and 0.4 respectively. Rapidity and transverse momentum cuts are identical to those 
used in panel (a) of these figures. 

We make the following observations on the basis of the two-particle A77 vs. A(j> correlations. First, we note that the 
near side "jet peak" tends to increase in amplitude progressively with larger radial velocity. The amplitude increase 
stems form the kincmatical focusing produced by the radial boost on the trigger particles and its associates: with 
large /3, more particles are focused on the near side peak in the given pt range of interest. 

Next, we observe the away-side ridge flattens out in azimuth at small velocity, progressively separates into two 
parallel structures with increasing velocity, and is refocused on the near side at the highest velocities considered here. 
We note the weak away-side "two-bumps" structures seen at low beta in projections on the axis are qualitatively 
similar to the structures reported by STAR and PHENIX based on two-particle correlations Q . The away-side ridge 
is essentially deflected and defocused away. , This leads to an increase of the correlation background for low f3r and to 
progressive build up of a ridge on the near side for larger boost velocities in qualitative agreement with observations 
by STAR |,[|. Radial flow can thus simultaneously produce the disappearance of the away-side "jet" peak, the 
formation of a dip structure, and the production of a strong ridge on the near side. 

We note in closing that given our simulations use one boosted p + p per simulated event, momentum conservation 
is globally violated. To achieve global momentum conservation, one would require addition of a recoiling background. 
Such background shall lead to a cos{A(j)) in the correlation functions thereby resulting in an increase of the away side 
yield [20j . The strength of this effect shall depend on event multiplicity, momentum ranges considered, and details of 
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FIG. 7: (Color online) Normalized two-particle correlations in A77 vs. Ai^ obtained with p+p PYTHIA events: Particle 1 (jet 
tag particle) in the range 3 < pt < 20 GeV/c, particles 2 and 3 (associates) in the ranges 2 < pt < 3 GeV/c. Panel (a) through 
(d) obtained with radial boosts /3r=0.0, 0.2, 0.3, 0.4 respectively. Middle and bottom subpanels show projections along A77 for 
\A(t>\ < 0.7 rad, and A(j> for |Ar;| < 1.0. See text for details. 

the collision dynamics. Effects of the recoil were therefore neglected in the simple model presented in this paper. 

IV. DISCUSSION 



The model used in this work is obviously simplistic. We do not expect parton-parton collisions in a dense-medium 
to proceed exactly as those in vacuum. We also do not expect all particles produced at a given radius in A-|-A to be 
subject to equal radial collective velocity. The use of radially boosted PYTHIA final states is thus not entirely justified. 
The spirit of the model should however be considered seriously as it provides a simple and elegant explanation to 
several observed phenomena simultaneously. A quantitative approach however requires attention to details: one must 
indeed consider modifications of the cross- section in the presence of the dense medium medium, gluon radiation, 
longitudinal flow effects, and rescattering or diffusion effects. 

Explaining the large energy and particle yield of the jet in terms of energy loss is also problematic given the 
observation that the ridge yield is basically independent of the selected trigger particle pt range, while the yield of the 
jet peak itself is independent of the selected pt. The radial boost mechanism, on the other hand, does not specifically 
modify the jet itself, but enables the formation of a strong ridge correlated in azimuth to the jet. The ridge thus 
carries high particle multiplicity yields, and energy, by construction. 
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FIG. 8: (Color online) Normalized three-particle azimuthal cumulant C3(A(^i2, Ai^ia) obtained with p+p PYTHIA events: 
Particle 1 (jet tag particle) in the range 3 < pt < 20 GeV/c, particles 2 and 3 (associates) in the range 2 < pt < 3 GeV/c. 
Panel (a) through (d) obtained with radial boosts /3r=0.0, 0.2, 0.3, 0.4 respectively. Middle and bottom panels show projections 
along cumulant diagonals. See text for details. 
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